Edge States and the Quantized Hall Effect in Graphene 
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We study edges states of graphene ribbons in the quantized Hall regime, and show that they can be 
described within a continuum model (the Dirac equation) when appropriate boundary conditions are 
adopted. The two simplest terminations, zigzag and armchair edges, are studied in detail. For zigzag 
edges, we find that the lowest Landau level states terminate in two types of edge states, dispersionless 
and current-carrying surface states. The latter involve components on different sublattices that may 
be separated by distances far greater than the magnetic length. For armchair edges, the boundary 
conditions are met by admixing states from different valleys, and we show that this leads to a single 
set of edges states for the lowest Landau level and two sets for all higher Landau levels. In both 
cases, the resulting Hall conductance step for the lowest Landau level is half that between higher 
Landau levels, as observed in experiment. 

PACS numbers: 73.43-f,73.20-r,73.23-b 



I. INTRODUCTION 

Recent progress in the processing of graphite have 
made possible the isolation of two dimensional carbon 
sheets, known as graphene^. This system has been stud- 
ied theoretically for a number of years, because when 
rolled up they form carbon nanotubes^. The material is 
unique because the underlying honeycomb lattice has a 
band structure with Dirac points at the corners of the 
Brillouin zone, two of which are inequivalent. Undoped, 
the system has one electron per atom and the Fermi en- 
ergy surface passes directly through the Dirac points. For 
low energies and dopings the system may be described by 
the Dirac equation. 

The stabilization of flat graphene sheets has allowed 
the application of perpendicular magnetic fields and the 
observation of the integer quantized Hall effect^. The 
striking result in these experiments is the first step height 
in the Hall conductance as a function filling factor, corre- 
sponding to filling the lowest Landau level (with electrons 
or holes), is half that of all subsequent steps. This be- 
havior was expected^ based on the bulk energy spectrum 
of graphene in a magnetic field: for a given spin species, 
there are pairs of Landau level bands at positive energies, 
each with partners at negative energies due to particle- 
hole symmetry. The lowest Landau level (LLL) has two 
levels precisely at zero energy, each of which is its own 
particle-hole conjugate^. This property of the LLL re- 
sults in its smaller contribution to the Hall conductance. 

The recent experiments on graphene studied ribbons 
that were relatively narrow, with widths in the microni 
or submicron- range. Under such circumstances, trans- 
port in the quantum Hall regime is typically dominated 
by edge states^. In this work we study edge states for 
graphene ribbons in detail, focusing on the simplest cases 
of a zigzag edge and an armchair edge. We demonstrate 
that a continuum description of edge states based on the 
Dirac equation is possible with the adoption of appropri- 
ate boundary conditions. In an edge state description. 



the quantization of the Hall conductance is determined 
by the number of edge state bands crossing the Fermi 
level. The Hall conductance results imply that the LLL 
supports only a single particle-like and a single hole-like 
band (per spin) at each edge, while the higher Landau 
levels have twice as many. Our goal is to understand how 
and why this happens, in a non-interacting picture. (Be- 
cause real spin plays no role in this study, we will from 
here on assume that all the electrons are spin-polarized, 
and refrain from explicitly noting the spin degree of free- 
dom in our discussions.) 

We now summarize our results. For zigzag edges, we 
show the correct boundary condition is for the wavefunc- 
tion to vanish on a single sublattice across the edge. In 
this case the LLL supports two types of edge states, 
which we call current-carrying and dispersionless sur- 
face states. Both states have strong components at the 
boundary of the system, but the former has equal weights 
on both types of sites of the honeycomb lattice, whereas 
the latter exists essentially only on one sublattice and has 
precisely zero energy. Such zero energy surface states are 
well-known to exist in graphene ribbons in the absence of 
a magnetic fieldi^, and have been shown in tight-binding 
calculations to persist when a field is applied^. In the 
quantum Hall context we find that the dispersionless sur- 
face states play a special role in forming two branches of 
edge states that do not pass through the Fermi level for 
any non-zero doping. The current-carrying edge states 
are also remarkable in that the surface contribution on 
one sublattice can be highly separated from a component 
well-inside the bulk of the sample, on the other sublattice. 
An interesting consequence of this is that a particle in- 
jected near a zigzag edge should oscillate back and forth 
between the edge and bulk, although presumably such 
oscillations would be damped by many-body effects not 
included in our study. 

For armchair edges, we find the correct boundary con- 
dition is vanishing of the wavefunction on both sublat- 
tices at the edge. This is achieved by mixing of wave- 
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functions from both Dirac points. In this case, there are 
no dispersionless surface states, and the LLL edge states 
behave differently than the higher Landau levels for other 
reasons. As we shall see, the energetics of states from one 
of these valleys is generically higher than from the other 
on a given sublattice, so that in the LLL only one band 
of edge states can meet the boundary condition, whereas 
in higher Landau levels there are two such bands. The 
admixing of the two valleys leads to wavefunctions with 
a characteristic spatial oscillation of the electron density 
with wavevector equal to AX^,, the difference between 
the Dirac point wavevectors along the direction perpen- 
dicular to the edge. Such oscillations should be observ- 
able in STM measurements. 



II. PRELIMINARIES 

We begin by reviewing some generalities about elec- 
trons in graphene. The lattice structure is a triangular 
lattice whose primitive lattice vectors are a = ao(l,0) 
and b — ao(l/2, \/3/2). There are two atoms per unit 
cell located at (0,0) and at d = ao(0, l/VS). A simple 
tight-binding model with only nearest neighbor hopping 
t leads to a Hamiltonian with Dirac points at the six 
corners of the Brillouin zone, only two of which are in- 
equivalent, and we take these to be K = a-nd 

~ ao'(~5' 'Tf-'' Wavefunctions can be expressed via 
the k-P approximation^iifi in terms of envelope functions 
[■(/'A(r), ■(/'-B(r)] and [t/;^ (r) , t/i^ (r)] for states near the K 
and K' points, respectively. These can be conveniently 
combined into a 4- vector — {^|JA, ips, (The 
reason for this sign convention will become apparent 
when we discuss the armchair edge.) This satisfies a 
Dirac equation H'i' = e^*, 



1/2, for /z = A,B. The solutions for states well away 
from the edge are well-known^. Taking A = —Bxy, 
^ = -^e^^yy^^ and V^(') = — ^e*'=''?'(^^(') with Ly 

the y extension of the sample, the wavefunctions retain 
their valley index as a good quantum number, and the 
positive energy wavefunctions may be written as $ = 
[(l)n-i{x-{ky+Ky)e),^r.{x~{ky+Ky)e),Q,Q\ior theK 
valley, and $ = [(),0, (t)n{x-{ky^K'y)l'^),-(j)n-i{x-(ky + 
K'y)^)] for the K' valley, with energies £„ = ^v^. In 
these expressions, (pn is the nth harmonic oscillator state 
and £ = ^JeBjc is the magnetic length. The negative 
energy states are easily obtained by reversing the signs 
of the wavefunctions on the B sublattice. For the case of 
the LLL, n = 0, and only one component of the 4- vectors 
for each valley is non-zero. This means the particle-hole 
conjugate of these wavefunctions are themselves. The 
bulk LLL wavefunctions do not have a clear particle- or 
hole-like character. 



H = 7ao 
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with 7 = \/3t/2. Note that 
in reciprocal space of the wavefunction from the K (K') 
point in the upper left (lower right) block of the Hamil- 
tonian. To apply this Hamiltonian in the presence of 
a magnetic field, one makes^ the Peierls substitution 
k — + — iV -I- eA/c where A is the vector potential. 

Before applying this procedure to systems with an 
edge, we point out some interesting and useful prop- 
erties of H. Firstly, H (and the more exact tight- 
binding Hamiltonian from which it descended) haa^ chi- 
ral (i.e., particle-hole) symmetry, THT — —H, where 
a, 
<Jz 

us that a solution to the Dirac equation ^E" with energy e 
has a particle-hole conjugate partner F^* with energy — e. 
Because of this, the wavefunctions must be normalized 
on each sublattice separately: / (ir[|-0^t(r)p -I- |V'u(i")P] = 



III. ZIGZAG EDGE 

The geometry for a zigzag edge is illustrated on the 
top and bottom edges of Fig. 1. It is interesting to note 
that each atom at the edge is of the same sublattice (say 
A). We shall see below that the appropriate boundary 
condition is to set the wavefunction to zero on a single 
sublattice (B), which we can understand to be the line of 
lattice sites that would lie just above or below the system 
if the bonds had not been cut to form the edge. In our 
discussions we will work with edges that lie along the y 
direction, so in what follows the coordinate axes in Fig. 1 
will be rotated by 90°. We begin by computing the band 
structure for a tight binding model of a graphene ribbon 
with zigzag edges, an example of which is illustrated in 
iky \ Fig. 2. The flat degenerate bands over a range of ky are 

I Landau levels, and, in the case of the LLL, dispersionless 

kx + iky I ' surface states which we discuss below. In a wide sam- 
kx — iky / pie, there is generically a large degeneracy within each 
(1) Landau band, because for the K (K') valley there are 
k denotes the separation wavefunctions peaked at Xp = [ky + Ky{') + nGy\f, 



and a-, is the Pauli matrix. This tells 



where Gy is a reciprocal lattice vector for the ribbon, 
and the integer n can take any value such that Xp is 
between the sample edges. One may conveniently reor- 
ganize the states by allowing all values of ky such that 
~Lx/2 < [ky + Ky{')]£'^ < Lx/2, with Lx the ribbon 
width, and assigning one state for each ky in the extended 
zone. 

The prominent structure in Fig. 2 is the appearance 
of dispersing energy bands, which occur when the wave- 
functions approach the edges. For the higher Landau 
levels one observes two pairs of such bands, whereas for 
the LLL there is only one such pair. This means that a 
Fermi energy crossing between the nth and (n-l-l)th Lan- 
dau levels yields a Hall conductance axy — (2n + l)e^//i, 
as observed in experiment2s4. 

The unique behavior of the LLL edge states may be 
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understood in terms of eigenstates of the Dirac Hamilto- 
nian with vanishing boundary conditions on a single sub- 
lattice. We begin by rotating the wavevectors in Eq. ^ 
kx ky, ky ^ —kx so that the zigzag edge lies along the 
y, and our wavefunctions then exist in the space a; > 0. 
Taking A = ~Bxy in this coordinate system, and defin- 



ing the ladder operator a 



-ky + x/£ + dx], with 



ky ky 

obey 



' KyC) for the K (K') valley, the wavefunctions 
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It is easy to see if one solves the equations for ip b and , 
the remaining wavefunctions are determined by (pA — 
a^ps/^ and f'g — —aip'j^/e. 

For the zigzag edge, the boundary condition does not 
admix valleys, and we can meet it for each type of wave- 
function separately: lpb{x — Q) = p'gix = 0) = 0. Thus 
for the K valley the spectrum e^{ky) is identical to that 
of a quantum Hall edge with a sharp boundarjiii. The 
wavefunctions ifB,n similarly are identical to their stan- 
dard Hall edge counterparts, turning into states in the 
7ith Landau level as kyi'^ moves well away from the edge. 
For n > 1, ipA,n is quite similar to a state in the (n — 1) 
Landau level provided the center of the wavefunction is 
not too close to the edge. 

In the LLL {n — 0), pA is qualitatively different. Be- 
cause LLL states are annihilated by the ladder operator 
a, and ipB is similar to a bulk LLL state when the cen- 
ter of the wavefunction is not too close to the edge, ipA 
is extremely small, except close to a; = where ips van- 
ishes and is forced to deviate from a bulk LLL state. The 
result is that (pA is strongly confined to the surface, and 
as the center of the tpB moves further into the interior 
of the sample, ipA becomes increasingly so confined. De- 
spite this strongly localized form, the normalization of 
the wavefunctions discussed in the previous section re- 
quires that fully half the probability of finding the elec- 
tron resides in this surface contribution. We note that, 
because e disperses with ky, these surface states carry 
current and contribute to the Hall conductivity. 

Within the Dirac equation, the existence and form of 
these current-carrying surface states can be examined 
with a variational approach. We adopt a trial wavefunc- 
tion (Pb{x) = w{x)(j)Q{x — kyP) and require w{x — > 0) = 

0. One may easily confirm that = 2 dx (^)'^ \(f)o\'^ , 
where e = {£/^/2jaQ)e. A simple choice for w is w = 



cs(l 



-Xx 



), where is a normalization constant and 



our variational parameter is A. With this choice one finds 
ifA^x) oc exp{— [a; — (ky — X)i'^]'^ /2P}, so that if A > ky, 
ipA is confined to the surface. The result of minimizing 
is illustrated in Fig. 3, where in the inset one sees that 
A increases faster than ky , so that tpA becomes more con- 
fined to the surface of the sample as ky increases and (pB 
penetrates into the bulk. 



Direct examination of wavefunctions from the tight- 
binding model confirms this basic picture, except in one 
important respect. Whereas the Dirac equation allows 
current carrying states with ipA increasingly localized to 
the surface as ky grows (and the peak position Xp of pB 
moves into the bulk) , the tight-binding results show that 
for Xp > \\Ky — K'y\P' the LLL wavefunction on sub- 
lattice A and the surface state on sublattice B appear 
as separate states, and that the surface state now moves 
back into the interior of the system with further increase 
in ky. The current-carrying surface state evolves into a 
dispersionless surface state which we describe below, and 
the number of allowed current-carrying surface states is 
limited. This is clearly an effect of the discreteness of 
the lattice that is not captured by the Dirac equation: 
independent, highly localized surface states can be writ- 
ten down for any ky in the continuum, but on the lat- 
tice states at ky and ky -I- Gy have the same periodicity 
along the y direction, and so cannot be made indepen- 
dent. Thus the number of current-carrying surface states 
is limited to ~ i^y\Ky — Ky\/A'K. For any practical pur- 
pose, this may be imposed by introducing a cutoff kc in 
the allowed fc^'s for the current-carrying surface states. 
For ky > kc, one may take the K valley wavefunction to 
have its bulk form. 

The edge state wavefunctions of the K' valley on the 
A sublattice are analogously identical to the well-known 
onesii, although is shifted upward by a single unit 
due to the ordering of the operators in the last 
of Eqs. 13 For the LLL, the bulk states oc (t)[){x — Xp) , 
p'g = 0, with X'p = kyP - K'yP, exactly satisfy Eqs. El 
with zero energy. Remarkably, these states are unaffected 
by the edge. Moreover, because Ky ^ Ky, there are val- 
ues values of ky in the extended zone for which X'p < 0, 
and the state is confined to the surface. These are the 
dispersionless surface states: they do not contribute to 
the Hall conductivity. The LLL wavefunctions of the 
tight-binding results around the center of the bands in 
Fig. 2(a) behave precisely as the Dirac equation results 
suggest: they are strongly confined to the surface, and 
continuously evolve into bulk LLL states on the A sub- 
lattice as ky increases. The dispersionless surface states 
supported by the LLL are the reason that it carries only 
half the Hall conductivity of the higher Landau levels for 
the zigzag edgei^. 



IV. ARMCHAIR EDGE 

The armchair edge is illustrated as the left and right 
edges Fig. 1, and the corresponding bandstructure from 
a representative tight-binding calculation appears in Fig 
2(b). Here the edge runs along the y direction, and no 
rotation of the figure is needed to represent our calcu- 
lations. Unlike the zigzag edge, the Landau bands all 
have dispersing states in the same regions of ky , but the 
LLL has one pair each of hole-like and particle-like edge 
states, while all the higher Landau levels have two. 
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To understand this from the viewpoint of the Dirac 
equation, we need to impose appropriate boundary con- 
ditions. In Fig. 1 one may see that the termination 
consists of a Une of A-B dimers, so it is natural to have 
the wavefunction amplitude vanish on both sublattices 
at X = 0. To do this we must admix valleys, and require 
^b{x = 0) = ip'g{x = 0) and ipAix = 0) = ^'j^{x = 0). 
Using the Dirac equation, and the fact that Ky = Ky, 
this second condition imphes dxfBlx^o = —dx'p'Blx^Q- 
To understand the effect of this on the solutions, it is 
convenient to combine the v'b's into a single wavefunc- 
tion defined for — oo < a: < oo: = ipb{x)9{x) + 
lp'b{—x)9{—x), with 9{x) the step function. The bound- 
ary conditions then amount to i^ix) and its derivative 
being continuous at x = 0. (This was the reason for 
our choice of relative sign in the 4- vectors of Section II.) 
From Eqs. |3 and |21 it is easy to see that V' obeys a 
Schroedinger equation [—d1 + U{x)]iIj{x) = e^^p{x), with 

U{x) = ^[{\x\/e - kyf - l/e + {2/e)e{-x)\. For 
large kyt"^, this double well potential, illustrated in Fig. 
4, has low energy states associated with the left well at 
e'^«3/2, 5/2,..., while for the right well one has states 
at « 1/2, 3/2, 5/2,.... We thus see there will be 
hybridization leading to pairs of edge states for all the 
higher Landau levels, whereas for the LLL there will be 
just a single such state. 

The admixing of different valley states to meet the 
boundary condition means that the wavefunction will os- 
cillate with period 271 /\Kx — K'^\. The behavior can ex- 
plicitly be seen in Fig. 5, which illustrates a LLL edge 
state from the tight-binding calculation. The apparent 
oscillation has precisely the period one expects for the 
valley mixing we introduced in the Dirac equation to 



meet the boundary condition. Although the period of 
this oscillation is very short (3.69A), it is in principle 
observable by STM measurements because the samples 
can be open to their environment^, in contrast to GaAs 
systems. 



V. CONCLUSION 

In this paper we have studied the edge states of 
graphene ribbons with zigzag and armchair terminations. 
We found in both these cases that a continuum descrip- 
tion in the form of the Dirac equation captures most fea- 
tures of the states found in tight-binding calculations, 
provided the wavefunction vanishes at the termination of 
the sample. For zigzag edges, we found the boundary 
condition can be met by wavefunctions within a single 
valley, leading to two types of edge states in the lowest 
Landau level, current-carrying surface states and disper- 
sionless surface states. The latter of these explains why 
the contribution to the quantized Hall coefficient from 
the LLL is only half that of higher Landau levels. For 
the armchair edge, we found that admixing of valleys is 
necessary to satisfy the boundary condition. For higher 
Landau levels, there are two hybridizations of the valley 
states for which this is possible, whereas in the LLL there 
is only one. This again leads to a contribution to the Hall 
coefficient from the LLL half the size of those from other 
occupied Landau levels. 
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nearest neighbor hopping) can introduce an energy depen- 
dence in the dispersionless surface states. However, the re- 
sulting edge states are not chiral and do not contribute to 
the Hall conductivity. See Ref. |3- 




FIG. 1: Ilustration of a graphene system with edges. Top and 
bottom edges are zigzag edges, left and right are armchair 
edges. 
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FIG. 2: Examples of energy bands for a graphene ribbon with 
periodic boundary conditions in the y direction and edges in 
the X direction. _B=100T (0.00126 flux quanta per unit cell.) 
Unit of energy ei = \^'yao/£. (a) Ribbon with zigzag edges, 
500 sites (530A) wide, (b) Ribbon with armchair edges, 1000 
sites (460A) wide. 
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FIG. 3: Surface part of current-carrying edge state for zigzag 
edge from variational method described in text. Main figure: 
ifA for two different dioices of ky. Inset: \mi-n (solid line) vs. 
ky (dashed line), demonstrating that (^a becomes increasingly 
localized on the surface with increasing ky. 



|U(x) 




FIG. 4: Potential U{x) for an armchair edge. See text. 
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FIG. 5: Squared wavefunction for an edge state of the arm- 
chair edge from tight binding calculation, ky = 2.l7r/ao and 
e/si = —0.202. Wavefunction penetrates sample over length 
scale I, while oscillations due to valley mixing occur on a much 
smaller length scale (3ao/2). 



